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ABSTRACT 


Context. Features in the spectra of primary cosmic rays (CRs) provide invaluable information on the propagation of these particles in 
the Galaxy. In the rigidity region around a few hundred GV, such features have been measured in the proton and helium spectra by the 
PAMELA experiment and later confirmed with a higher significance by AMS-02. We investigate the implications of these datasets 
for the scenario in which CRs propagate under the action of self-generated waves. 

Aims. We show that the recent data on the spectrum of protons and helium nuclei as collected with AMS-02 and Voyager are in very 
good agreement with the predictions of a model in which the transport of Galactic CRs is regulated by self-generated waves. We 
also study the implications of the scenario for the horon-to-carhon ratio: although a good overall agreement is found, at high energy 
we find marginal support for a (quasi) energy independent contribution to the grammage, that we argue may come from the sources 
themselves. 

Methods. The transport equation for both primary and secondary nuclei is solved together with an equation for the evolution of the 
self-generated waves and a background of pre-existing waves. The solution of this system of non-linear equations is found with an 
iterative method elaborated by the same authors in previous work on this topic. 

Results. A break in the spectra of all nuclei is found at rigidity of a few hundred GV, as a result of a transition from self-generated 
waves to pre-existing waves with a Kolmogorov power spectrum. Neither the slope of the diffusion coefficient, nor its normalisation 
are free parameters. Moreover, at rigidities below a few GV, CRs are predicted to be advected with the self-generated waves at the 
local Alfven speed. This effect, predicted in our previous work, provides an excellent fit to the Voyager data on the proton and helium 
spectra at low energies, providing additional support to the model. 


1. Introduction 


The transport of Galactic cosmic rays (CRs) is likely to be very 
complex: the structure of the large scale magnetic field is compli¬ 
cated and very poorly known, and the structures on small scales 
(resonant with the CR particles) are basically unknown, although 
the fact that the diffusion paradigm seems to work can be con¬ 
sidered as an indirect evidence for the existence of such small 
scale turbulence. The origin of the power on small scales is also 
unknown, but one contribution that seems to be hardly avoid¬ 
able is that due to the self-generation of perturbations due to the 
CR current, proportional to the gradient of CRs, which in turn 
is due to the existence of the same scattering centres, responsi¬ 
ble for diffusion. This simple description is sufficient to empha¬ 
sise the non-linear nature of this process, which is qualitatively 
similar to what happens at supernova shocks, thought to be the 
main s ources of Galactic CRs (see the recent review by iBlasi I 
(l2Q13l) ). In addition to the self-generated waves, turbulence at 
larger spatial scales is generically expected and becomes impor¬ 
tant for scattering CRs of higher energies. 

This scenario has been analysed in detail in (iBlasi et al. I 


1201 21: lAloisiol 20131), whe re the main implications were dis- 
cussed: IBlasi et al!T(l2012h calculated the spectrum of protons 


under the action of both the self-generated and pre-existing tur¬ 
bulence, and co mpared the res ults with the PAMELA data avail¬ 
able at the time dAdrian il201 Ih . where the first direct detection of 
a spectral break in proton and helium fluxes at few hundred GV 
was claimed. The first release of the data collected by AMS-02 
(iHaino 11201 3t IChoutko1l2013l) did not confirm the existence of 
these spectral features and brought the investigation on this topic 
to an almost complete standstill, waiting for the resolution of 
the observational conundrum. Recently the AMS collaboration 
published the final a nalysis of the data on the proton spectrum 
(lAguilar et al.ll2015h . where a change of slope at few hundred 
GV is evident. At the present time, only preliminary results on 
the sp ectrum of helium and carbon nuclei are available (lAMS-021 
l2015l) . but a similar break is visible in such data as well. More¬ 
over, preliminary data on the B/C ratio have also been presented: 
the small statistical error bars up to high energies allow us to 
use this tool as a powerful indicator of the propagation of CRs 
through the Galaxy. 

In addition to the AMS-02 data, the results of another invalu¬ 
able experimental effort became available in the last few years: 
the Voyager spacecraft, launched in 1977, reached the termina¬ 
tion shock and is now believed to be mov ing in the interstel - 
lar medium, unaffected by the solar wind (iStone et al. II2Q13I) . 
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The Voyager is therefore providing us with the very first mea¬ 
surement of the interstellar spectra of protons and helium nuclei, 
that can be compared directly with the predictions of our models, 
rather than dealing with complex and uncertain recipes of solar 
modulation. 

These recent developments revived the interest of the Com¬ 
munity and stimulated a new search for explanations of the spec¬ 
tral brea ks: in this paper w e reconsider the model first put for¬ 
ward bv iBlasi et al. I (l2012h : lAloisiol (1201 3h and check it versus 
both the data of AMS-02 and Voyager. We find that the set of 
parameters that were previously used to fit the PAMELA data 
also lead to a good fit to the AMS-02 data on the proton and 
helium spectra. The rigidity where the spectral break occurs is 
predicted rather solidly by the model. In fact the model predicts 
in a rather constrained manner both the normalisation and the en¬ 
ergy dependence of the diffusion coefficient at rigidities below a 
few hundred GV, where self-generation is important. At the same 
time it also fixes the injection rate (for instance in SNR shocks) 
and hence the diffusion coefficient at high energies, in order to 
reproduce the observed spectrum of CRs. Since both contribu¬ 
tions are quite well constrained, so is the transition energy. The 
constrained nature of the model is also its main strength. As 
discussed above, the B/C ratio was prese nted by the AM S-02 
collaboration only in a preliminary form (lAMS-02ll2015h and 
we use this information only to further check, at least visually, 
the viability of the model: the AMS data extend to unprecedent¬ 
edly high energies and hence they provide us with much infor¬ 
mation about propagation. We find that the preliminary data on 
B/C are well fit with our model, although there is marginal ev¬ 
idence for some additional grammage with a harder energy de¬ 
pendence compared with the grammage accumulated through¬ 
out the Galaxy. The necessary grammage is compatible with that 
accumulated by CRs while being accelerated and advected in a 
typical SNR, and is consistent with being energy independent. 

The paper is structured as follows: in ^we summarize th e 
essential aspec ts of the model put forward by Blasi et all d2012l) : 

oisi ^ (I2013h . In 0we discuss the implications for the nuclear 
spectra and for the B/C ratio. We summarize in ^ 


2. Self-Generated versus pre-existing waves and 
CR transport 

In a ID model (infinite slab, only dependence from distance z 
above/below the plane retained), the transport equation for CR 
nuclei in its general form can be written as follows: 


stressed above, for the sake of a meaningful comparison with 
data, it is important to take into account the stable isotopes of 
all elements. This is important not only for pure secondary ele¬ 
ments, namely elements produced only through spallation, such 
as Boron (B='°B-i-^'B), but also for Nitrogen (N^^'^N-i-'^N), 
which gets a significant secondary contribution from spalla¬ 
tion of heavier nuc lei, for Carbon (C=*^C - H*^C) and Oxyge n 
(0=*®0 -h1^0h-'^ 0) (Idi Bernardo et al. ll^^lEvoli et al. II 200 I . 
Moreover, for simplicity we assume an instantaneous decay for 
isotopes whose lifetime is much shorter than their escape time 
from the Galaxy. This means that in the sum over a' in the rhs of 
Eq. ([1]) we consider also terms of the type cra'^a"fa' being a” a 
nuclear specie that rapidly decays into a. 

Here v{p) - ffc 'k the velocity of nuclei of type a having 
momentum p. Notice that since the gas is assumed to be present 
only in the disc, and the ionization rate is proportional to the 
gas density, one can write: (^) = 2hd6{z)bQdp), where 

bo.rtip) contains the particle p hysics aspects of the process (see 
IStrong & Moskalenko I (Il998l) and references therein for a more 
detailed discussion of this term). 

Eollowing the procedure outlined bv lAloisi^(l2013h . one can 
transform Eq. dB in a modified weighted slab transport equation: 


IJE) ^ d 
XdE) dE 



crglgiE) 

m 


AgP^qO^g(p) ^ Ig(E) 

= 2hd --H > ,-cr^'^Q., (2) 

u V m 

' a’>a 

where /„(£’) is the flux of nuclei with kinetic energy per nucleon 
E for nuclei of type a, such that Ia{E)dE = v p^fo^a{p)dp. It is 
easy to show that Ia{E) = Aap^fo^a(p), being Ag the atomic mass 
number of the nucleus. We also introduced the quantity 


Xg{E) = ^ 
2va 




(3) 


that represents the grammage for nuclei of type a with kinetic 
energy per nucleon E, while 


dE\ 

d^lad 


2va 




(4) 


is the rate of adiabatic energy losses due to advection. 

The diffusion coefficient relevant for a nucleus a can be writ¬ 
ten as: 
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( 1 ) 


where (jg is the spallation cross section of a nucleus of type a, 
// is a grammage parameter fixed to 2.4 mg/cm^, and qo^aip) 
is the rate of injection per unit volume in the disc of the 
Galaxy. Namely, since hd is the half-thickness of the (assumed 
infinitesimal) gaseous disk, 2hdqo^a is the rate of injection in 
the disc of the Galaxy per unit surface. The total cross section 
for spallation and the cross sections for the individual chan¬ 
nels of spallation of a heavi er element to a lighter elem ent 
{(Tg’-,g) have been taken from I Webber et alTI d 199(1 1^031) . As 


where W(k) is the power spectrum of waves at the resonant 
wavenumber k - ZgeB^fpc. The non-linearity of the prob¬ 
lem is evident here: the diffusion coefficient for each nuclear 
species depends on all other nuclei through the wave power 
W(k), but the spectra are in turn determined by the relevant dif¬ 
fusion coefficient. The problem can be written in a closed form, 
though in an implicit way, by using the transport equation for 
each n ucleus, Eq. dB, and writi ng the evolution equation for the 
waves (iMiller & R^erts]ll995ll : 


d dW 

- Ejik - 

dk dk 


+ TcrW - qw(k). 


(6) 


where qw(k) is the injection term of waves with wavenumber 
k. In the present calculations we assume that waves are only 
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injected on a scale Ic ~ 50 - 100 pc, for instance by super¬ 
nova explosions. This means that qwik) oc 6(k - I He)- The 
level of pre-existing turbulence is normalized to the total power 
rjB — SB^IBq - J dkW(k). Strictly speaking the wave number 
that appears in this formalism is the one in the direction parallel 
to that of the ordered magnetic field. In a more realistic situation 
in which most power is on large spatial scales, the role of the 
ordered field is probably played by the local magnetic field on 
the largest scale. 

The term TcrIT in Eq. ® describes the generation of wave 
pow er through CR induced streaming instability, with a growth 
rate (ISkilling Ill975h : 


rcr(^) = 


167r^ 


VA 


k W(k)Bl 


z 


4 . .dfa 
P 


(7) 


- p=ZaeBQlkc 


where a is the index labeling nuclei of different types. All nu¬ 
clei, including all stable isotopes for a given value of charge, 
are included in the calculations. As discussed in much previous 
literature, this is very important to compute properly the diffu¬ 
sion coefficient and thus for a meaningful comparison with the 
flux spectra and secondary to primary ratios, notably B/C. The 
growth rate, written as in Eq. O, refers to waves with wave num¬ 
ber k along the ordered magnetic field. It is basically impossible 
to generalize the growth rate to a more realistic field geometry 
by operating in the context of quasi-linear theory, therefore we 
will use here this expression but keeping in mind its limitations. 

The solution of Eq. (|6]l can be written in an implicit form 


W(k) = 





H- 


1 + a'2 
CkVa 


r” dk' 

. . . 1 

L ^ 

dkrcR(k)w(k) 

Jko 


(8) 


being ko - 1/4. In the present paper we assume a Kolmogorov 
phenomenology for the cascading turbulence, so that ai - 1 jl 
and a 2 - 1/2, and an unperturbed magnetic field Bq = l^G. 
The two terms in Eq. ([8]l refer respectively to the pre-existing 
magnetic turbulence and the CR induced turbulence. In the limit 
in which there are no CRs (or CRs do not play an appreciable 
role) one finds the standard Kolmogorov wave spectrum 


W(k) = Wo 



ai - I 
Q'2 + 1 


5 

3 


(9) 


normalized, as discussed above, to the total power VTo = (■? - 

1)4?7b- 

The equations for the waves and for CR transport are solved 
together in an iterative way, so as to return the spectra of par¬ 
ticles and the diffusion coefficient for each nuclear species and 
the associated grammage. The procedure is started by choosing 
guess injection factors for each type of nuclei, and a guess for 
the diffusion coefficient, which is assumed to coincide with the 
one predicted by quasi-linear theory in the presence of a back¬ 
ground turbulence. The first iteration returns the spectra of each 
nuclear specie and a spectrum of waves, that can be used now 
to calculate the diffusion coefficient self-consistently. The pro¬ 
cedure is repeated until convergence, which is typically reached 
in a few steps, and the resulting fluxes and ratios are compared 
with available data. This allows us to renormalize the injection 
rates and restart the whole procedure, which is repeated until a 
satisfactory fit is achieved. Since the fluxes of individual nuclei 
affect the grammage through the rate of excitation of stream¬ 
ing instability and viceversa the grammage affects the fluxes, the 
procedure is all but trivial. 



Fig. 1. Spectrum of protons measure d by Voyager (blue empty cir¬ 
cles), AMS-02 (black filled circles) dAguilar et al]|2015l) , PAMELA 
(green empty squa res) dAdrianil2011h and CREAM (blue filled squares) 
dYoon et al. Il201 ih . compared with the prediction of our calculations 
(lines). The solid line is the flux at the Earth after the correction due to 
solar modulation, while the dashed line is the spectrum in the ISM. 



Fig. 2. Spectrum of He nuclei according to preliminary measurements 
of AMS-02 (black filled circles), as measured by Voyager (blue empty 
circles), PAMELA (green empty squares) and CREAM (blue filled 
squares), compared with the prediction of our calculations (lines). The 
solid line is the flux at the Earth after the correction due to solar modu¬ 
lation, while the dashed line is the spectrum in the ISM. 


3. Results 

The main evidence for a transition from self-generated waves to 
pre-existing turbulence can be searched for in the spectra of the 
light elements, protons and helium nuclei. A spectral br eak was 
in fact found by the PAMELA experiment (lAdrianil201 Ih in both 
spectra and later confirmed by AMS-02, although at the time 
of writing this paper only the resu lts of AMS on protons have 
been published (I Aguilar et al.l2015l) . while a pre liminary versio n 
of the spectrum of helium has been presented (lAMS-02ll2015ll . 
The spectra of both elements were also measured by the Voyager 
(IStone et al. II2013I) outside the heliosphere, so as to make this 
the first measurement in human history of the CR spectra in the 
interstellar medium. This is a very important results in that it 
also allows us to refine our u nderstanding of the effects of solar 
modulation (lPotgieterll2Q13h . 

The spectrum of protons and helium nuclei as calculated 
in this paper is shown in Pigs. [1] and |2] respectively: the solid 
lines indicate the spectra at the Earth, namely after s olar mod- 
ulation modelled using the force-free approximation dGleeson I 
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Il968h . while the dashed lines are the spectra in the ISM. The 
data points are the spect ra measured by the Voyag er (empty cir¬ 
cles) l^one et al. Il2013h . AMS-02 (filled circles) (lAguilar et al.l 
l2015h . PAMEL A (empty squares) (lAdrianill201 ih and CREAM 
(filled squares) (lYoon et al. 11201 Ih . Eigs. [1] and |2] show several 
interesting aspects: 1) both the spectra of protons and helium nu¬ 
clei show a pronounced change of slope at few hundred GeV/n, 
where self-generation of waves becomes less important than pre¬ 
existing turbulence (in fact, the change of slope takes place in 
rigidity). 2) We confirm that injecting He with a slightly harder 
spectrum with respect to protons (p^”* versus improves 

the fit to the data. 3) The spectra calculated to optimise the fit to 
the AMS-02 and PAMELA data is in excellent agreement with 
the Voyager data (see dashed lines). This is all but trivial: in our 
model, at sufficiently low energies (below ~ 10 GeV/n), particle 
transport is dominated by advection (at the Alfen speed) with 
self-generated waves rather than diffusion. This reflects into a 
weak energy dependence of the propa gated spectra tha t is ex¬ 
actly what Voyager measured (see also (IPotgieter ll^oT^ f. 4) At 
low energies, the agreement of the predicted spectra with those 
measured by Voyager is actually better than the agreement with 
the modulated spectra, as observed with AMS-02; this suggests 
that probably the prescriptions used to describe solar modula¬ 
tion are somewhat oversimplified, either when applied to data 
collected over extended periods of time, when the effective solar 
potential may change appreciably, or because of intrinsic limita¬ 
tions of the force-field approximation. 

Eor each heavier nucleus, we assume the same injected spec¬ 
tral shape in rigidity as for helium, keeping as only free param¬ 
eter the normalization, chosen to match the data. In Eig.|3]we il¬ 
lustrate the prediction for Carbon nuclei (which is also a needed 
ingredient to compute the B/C ratio), compared with data by 
PAMELA and CREAM, as well as preliminary data by AMS- 
02. The free normalization is chosen to match more closely the 
AMS-02 data. Clearly, the phenomenon of transition from self¬ 
generated to pre-existing waves manifests itself in the transport 
of all nuclei, hence we should expect a spectral break at the 
same rigidity as for helium and protons. This prediction appears 
currently in agreement with Carbon spectrum observations, al¬ 
though it is hard to judge to what extent a break is present in 
AMS-02 data alone, giving the growing error bars and the lim¬ 
ited dynamical range at high energy. Note that a break would 
appear more prominent if one were to combine PAMELA and 
CREAM data, which do seem to differ from AMS-02 data in the 
10 to ~ 200 GeV/n range beyond the reported errors. Definitely, 
the forthcoming AMS-02 publication of nuclear fluxes should 
help in clarifying the situation. 

In Eig.|4]we show the calculated B/C ratio (solid black line) 
as compared with data from CREAM (blue squares), PAMELA 
(green squares), and the still preliminary ones from AMS-02 
(black circles). Even if the injected Carbon flux is normalized 
to the preliminary Carbon data reported by AMS-02, the B/C 
ratio is still in satisfactory agreement with both PAMELA and 
CREAM data, as for our previous result (lAloisiol l2013h . The 
B/C ratio also fits the AMS data up to ~ lOOGeV/n. At higher 
energy, the AMS-02 analysis seems to suggest a B/C ratio some¬ 
what higher than our prediction. While its significance is uncer¬ 
tain, given the preliminary nature of AMS data, if this “excess” 
is interpreted as physical, it would suggest the presence of an 
additional contribution to the grammage traversed by CRs. The 
most straightforward possibility to account for such a grammage 
is that it may be due to the matter traversed by CRs while es¬ 
caping the source, for instance a SNR. The grammage due to 



Fig. 3. Spectrum of C nuclei as measured by CREAM (blue squares), 
PAMELA (green empty squares), and according to preliminary mea¬ 
surements of AMS-02 (black circles), compared with the prediction of 
our calculations (lines). The solid line is the flux at the Earth after the 
correction due to solar modulation, while the dashed line is the spectrum 
in the ISM. 



Fig. 4. B/C ratio as measured by CREAM (blue squares), PAMELA 
(green empty squares), and according to preliminary measurements of 
AMS-02 (black circles). The black/bottom solid line is the prediction of 
our model, while the red/top line has been obtained by adding a source 
grammage of 0.15 g cm“^, close to that given by Eq. flOl l. 

confinement inside a SNR can be easily estimated as 

^SNR - 1.4r,m^/7isMcrsNR - 0.17 , (10) 

cm“-^ 2 X lO^yr 

where hism is the density of the interstellar gas upstream of a 
SNR shock and r, = 4 is the compression factor at the shock 
and Tsnr is the duration of the SNR event (or better, the lifetime 
“useful” to confine particles up to E ~TeV/n), assumed here to 
be of order twenty thousand years. The factor 1.4 in Eq. (fTOl) has 
been introduced to account for the presence of elements heavier 
than hydrogen in the target. While Eq. (fTOl) is only a rough es¬ 
timate of the grammage at the source, in that several (in general 
energy dependent) factors may affect such an estimate, at least 
it provides us with a reasonable benchmark value. The solid red 
curve in Eig. |4] shows the result of adding the grammage accu¬ 
mulated by CRs inside the source to the one due to propagation 
in the Galaxy. It is clear that by eye it fits better the AMS-02 
data at high rigidity, while being also compatible with the older 
CREAM data. The forthcoming publication by AMS-02 of the 
fluxes of nuclei and secondary to primary ratios should hope- 
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fully clarify the situation and provide more reliable discriminat¬ 
ing power. 

4. Summary and Conclusions 

One of the main implications of the transport theory of CRs in 
the Galaxy is that the spatial gradient developed due to the dif¬ 
fusive propagation of CRs (for instance tuned to fit the B/C ra¬ 
tio) is such that the rate of growth of Alfven waves cannot be 
neglected, and in fact at rigidity < 200 GV the self-generated 
waves must play a dominant role. On the other hand, at higher 
rigidity the growth rate drops with respect to damping processes 
responsibl e for cascad i ng of waves in wav enumber space . As 
shown by iBlasi et ahl d2012ll and later by I Alois! 3 (1201 3ll the 
transition between the two reflects in a change of slope in the 
spectra of primary elements, most notably protons and helium 
nuclei. Evide nce in favou r of such a break emerged first in nu¬ 
clear spectra (lAhnll201^ : then hints for something similar in 
proton and helium spectra were suggested by the need to recon- 
cile high-energy data measured by CREAM (I Yoon et al. 1 1201 1) 
with low-energy ones; more recently, PAMELA (lAdrianIbol 1 ) 
reported for the first time the detection of the transition in pro¬ 
ton and helium nuclei; Anally, after contradictory early result s 
presented by AMS-02 in 2013 (iHaino 1 120 1 31: IChoutko 1 120 1 3l) . 
this behaviour for proto n and helium has been confirmed by the 
AMS-02 collaboration (I Aguilar et al.ll20L5l: IAMS-02|IM1 . al¬ 
beit the publication only includes the proton spectrum, with the 
helium flux being still preliminary. 

The results found in this paper can be summarised as fol¬ 
lows: 1) both the proton and helium spectra show a pronounced 
change of slope at rigidity of a few hundred GV, which we inter¬ 
pret as the region where the self-generation of waves becomes 
less important than pre-existing turbulence. Our interpretation 
implies that a similar spectral break should also be present in 
the spectra of other primary nuclei (for instance carbon); cur¬ 
rent data are consistent with this prediction, although the un¬ 
certainties, both the ones of statistical nature and possibly oth¬ 
ers of systematic origin, do not enable a stringent quantitative 
test, so that no final conclusion on this issue can be drawn, 
yet. 2) The observed difference in spectral slope between He 
nuclei and protons should be attributed to differences at accel¬ 
eration/injection, with the injection spectrum of He nuclei that 
must be somewhat harder (cx p '^ '^) compared with the pro¬ 
ton one (oc p "^ ^). A consensus on the theoretical explanation of 
this difference is not yet available, although some models have 
been put forward, for instance involving: different a cceleration 
sites (as in the reverse shock of iPtuskin et al. I (l2013h . the inho¬ 
mogeneous chemical compositio n of the medium at the accel¬ 
eration site (IOhira&lokall2ni1h . o r phenomena associa ted to 
the physics of the injection process (iMalkov et al. ll2()T^ . such 
as that the more easily ionized hydrogen gets a greater contri- 
buti on by older a nd weaker shocks (associated to steeper spec¬ 
tra) dPrurv 11201 Ih . 3) The spectra of protons and He nuclei pre¬ 
dicted by our model automatically fit the low energy Voyager 
data, which is remarkable, because this suggests that the propa¬ 
gation of low rigidity (< 10 GV) protons and He nuclei is dom¬ 
inated by advection wi th self-generated A lfven waves, as was 
originally proposed by iBlasi et ^ (l2012h . 4) Data on secon¬ 
daries or secondary to primary ratios pre- AMS are consi stent 
with the picture above, as already shown in (lAloisioll2013h . but 
do not allow for a stringent test of this scenario, yet. Preliminary 
AMS-02 data on the B/C ratio show greater potential for discrim¬ 
inating power and, if taken at face value, provide marginal evi¬ 
dence for an excess grammage at energy/nucleon >100 GeV/n. 


We speculated that this grammage may be the one traversed by 
CRs in the acceleration region: in fact, a quick estimate of the 
equivalent grammage in a typical SNR (putative source of CRs) 
compares exceptionally well with the grammage needed to rec¬ 
oncile the predicted and the observed one. A generic expecta¬ 
tion is that at Energies/nucleon of 0(TeV/n), the contribution 
of production at the “sources” for so-called secondary species 
may become significant, as already noted in the p ast for other 
species such as antiprotons (IBlasi & Serpico ir2009t) . Needless to 
say, current B/C data presented by the AMS-02 collaboration 
are still preliminary and one should not over-interpret them, es¬ 
pecially in the light of the challenging path that lead to the Anal 
establishment of the presence of changes of slope in the proton 
(and most likely also helium) spectrum. Yet, it is certainly true 
that secondary-to-primary ratios, at the level of precision and 
dynamical range that should be achievable by AMS-02, carry a 
substantial amount of information, as illustrated by our compar¬ 
ison. A more detailed re-analysis of our model will be certainly 
justified, once AMS-02 data on nuclear species (both primaries 
and secondary ones) will be eventually published. 
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